Detectability of tensor modes in the presence of foregrounds 



Mihail Amarie,"'^'Q Christopher Hirata,^ and Uros Seljak^' ^ 

^Department of Physics, Jadwin Hall, Princeton University, Princeton, New Jersey 08544, USA 
^International Center for Theoretical Physics, Strada Costiera 11, 34OI4 Trieste, Italy 

(Dated: February 5, 2008) 

In inflationary models gravitational waves are produced in the early universe and generate B-type 
polarization in the cosmic microwave background (CMB). Since B polarization is only generated 
by gravity waves it does not suffer from the usual cosmic variance. A perfect decomposition of 
the CMB into B-modes and _E-modes would require data from the entire sky, which in practice 
is not possible because of the foreground contaminants. This leads to mixing of E polarization 
into B, which introduces cosmic variance contamination of B polarization and reduces sensitivity 
to gravity wave amplitude even in absence of detector noise. We present numerical results for the 
uncertainty in the tensor-to-scalar ratio using the Fisher matrix formalism for various resolutions 
and considering several cuts of the sky, using the foreground model based on dust maps and assuming 
90GHz operating frequency. We find that the usual scaling A (-^^ oc f~i}J^ is significantly degraded 

and becomes A (-^j oc /J^^ for fsky > 0.7. This dependence is affected only weakly by the choice 
of sky cuts. To put this into a context of what is required level of foreground cleaning, to achieve a 
T/S = 10~^ detection at 3 a one needs to observe 15% of the sky as opposed to naive expectation of 
0.3%. To prevent contamination over this large sky area at required level one must be able to remove 
polarized dust emission at or better than 0.1% of unpolarized intensity, assuming the cleanest part 
of the sky has been chosen. To achieve T/S — 10"'* detection at 3 ct one needs to observe 70% 
of the sky, which is only possible if dust emission is removed everywhere over this region at 0.01% 
level. Reaching T/S = 10"^ should be easier: 1% of the sky is needed over which polarized emission 
needs to be removed at 1% of unpolarized intensity if the cleanest region is chosen. These results 
suggest that foreground contamination may make it difficult to achieve levels below T/S = 10^^. 

PACS numbers: 98.70.Vc, 98.80.Bp, 98.80.Cq 



I. INTRODUCTION 



The polarization in the CMB is generated by Thomson 
scattering of the CMB photons off free electrons. Thom- 
son scattering generates only linear polarization and the 
polarization can be expressed in terms of the coordinate 
dependent Stokes parameters Q and U. It can be also 
decomposed in the coordinate independent components 
in the harmonic space: scalar components denoted by E, 
and pseudo-scalar components B. For spin I, the parity of 
the i?-modes is (—1)', while the parity of the i3-modes is 
— (— 1)'. The primordial density perturbations are scalar 
only, and because of the parity invariance, they can only 
generate E modes. The gravitational waves present dur- 
ing the inflation are tensor perturbations, so they can 
generate both E and B modes til i3| • If the ampli- 
tude of the gravity waves is very small relative to scalars 
it cannot be isolated from the temperature anisotropics 
or E polarization due to cosmic variance. The B po- 
larization is insensitive to the cosmic variance from the 
E modes, being affected only by the foregrounds and 
the instrument noise. The amplitude of the B modes 
depends on the amplitude of the gravity waves gener- 
ated during the inflation, which in turn depends on the 
(as yet unknown) energy scale at which the inflation oc- 
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curred. The tensor to scalar power ratio is often de- 
fined asT/S = /Q^'^'^'^^^'^'^ ^ i g_ ratio of the 
contributions to the CMB temperature quadrupole. In 
terms of the energy density during the inflation it can 
be written as T/S - V;/(3.7 x IQi^CeV)*, where K is 
the energy scale of inflation. The possibility of detect- 
ing direct evidence of gravity waves from inflation via 
_B-modes is currently being considered by a number of 
ground, balloon and space based experiments. Notably, 
a future satellite mission dedicated to B type polarization 
has been identifled as one of the NASA Beyond Einstein 
missions, the Inflation Probe. 

There are a number of outstanding issues that need to 
be resolved for such a mission to become a reality. Some 
of these are the question of required sensitivity and angu- 
lar resolution to achieve the desired energy scale of infla- 
tion below 10^^ GeV 3 and the issue of B-mode contam- 
ination from gravitational lensing j^. One of the main 
issues related to the feasibility of such observation of the 
CMB polarization is the contamination by foregrounds. 
The dust present in the interstellar medium (ISM) is one 
of the main foreground contaminants. It has a thermal 
vibrational emission spectrum, which is increasing with 
the frequency in the range of interest for us, and it also 
shows an increased emission in the 10-90 GHz range, 
which could be caused by the rotation of small dust par- 
ticles 0. The ISM can also be a significant source of 
synchrotron radiation, which shows a red frequency de- 
pendence 0- The ISM contaminants are present in our 
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galaxy, so the most affected region on a map would be 
the Galactic Plane. 

All these contaminants lead to a more realistic scenario 
in which we have to exclude portions of the sky because 
of their contamination by foregrounds. This can be done 
by simply removing some portions of the sky where the 
contamination is largest, or can be more sophisticated 
by marginalizing over some foreground templates 0, 0| . 
For simplicity, in this paper we only consider the first 
case. We expect that our conclusion hold qualitatively 
also for other more sophisticated types of analysis and 
discuss what are required conditions when this is likely 
to be valid. Most CMB parameter forecasting studies 
have assumed that if only a fraction of the sky /sky is ob- 
served then the uncertainties in all parameters is reduced 

— 1/2 

by /gjjy . However, in case of polarization this ignores 
an important issue. On the whole sky, the E and the 
B modes are separable. This statement does not hold 
anymore on a partial sky as the boundaries of the cuts 
generate mixing between these modes. In simple geome- 
tries one can isolate pure E modes, pure B modes and 
also mixed modes for which one cannot discern the E or 
B nature |lOl ITU . Pure E modes and mixed modes are 
contaminated by E mode polarization from scalar pertur- 
bations, which are larger than the B modes from gravity 
waves even at the present upper bound on T/S. Thus 
for these modes some of the same cosmic variance issues 
arise as for E polarization and temperature anisotropics. 
As a result the scaling with /sky can become much worse 
than expected. A few studies of tensor B-modes have 
taken this issue into account in specific cases such as the 
WMAP KpO cut or cuts based on Galactic latitude 
[T^ flSj . but the full dependence on the sky cut has not 
been investigated. 

Analytic decompositions into E, B and mixed modes, 
emphasized in previous work, are possible only in a few 
cases of simple geometry. In a more general case with pix- 
elization effects included, all modes are mixed, but the 
level of mixing varies so that the E mode contamination 
ranges from negligible to nearly complete. Qualitatively, 
the larger the area of the sky coverage and the smaller 
the length of the boundaries, the smaller the contamina- 
tion and mode mixing, but analytic results are difficult 
to obtain. Fortunately the decomposition into pure and 
mixed states, while useful for heuristic interpretation, is 
not necessary in the actual analysis of the data. One 
simply needs to perform the usual likelihood analysis of 
the data expressed with measured Stokes Q and U pa- 
rameters. The likelihood analysis is an optimal analysis 
and as such cannot be improved upon by performing the 
decomposition into pure and mixed modes. 

Broadly speaking there are two contributions to B- 
mode polarization. One is from the recombination epoch 
(which defines the last scattering surface), which leads 
to a peak in the power spectrum at I ~ 100. The sec- 
ond is from the reionization, which re-scatters the CMB 
quadrupole on the horizon scale during reionization and 
leads to a large signal at large angular scales, I < 20, 



while at the same time suppressing the recombination 
signal at smaller scales [l^. If the optical depth to reion- 
ization is high, as suggested by WMAP first year results 
E3iI3IS'la|) then the latter signal dominates in the de- 
termination of the gravity wave amplitude . Because of 
its larger angular scale, this component of the signal also 
suffers more from i?-mode mixing into _B-modes due to 
incomplete sky coverage. In this paper we include both, 
varying the amplitude of optical depth over the allowed 
range. 



II. QUADRATIC ESTIMATOR AND FISHER 
MATRIX 



We consider the quadratic estimator for the vector 
X = {Qi,Q2, ■■■Qn,Ui,U2, ■■■Un), where N is the total 
number of pixels considered and Q and U are the co- 
ordinate dependent Stokes parameters. The covariance 



matrix C = 



is written as 



T 



(1) 



Here Cq is the covariance matrix including noise and 
scalar fluctuations but no tensors, T/S is the tensor-to- 
scalar ratio, and Ct/s is the _B-mode power spectrum 
template. 



fir<BB 



(2) 



We have left out the -E-mode power spectrum from the 
tensors, which contributes negligibly at small T/S due 
to cosmic variance. In our case, the Co matrix contains 
both the contribution from the detector noise and the 
i?-mode contribution to the polarization. We work in 
the limit of T/S = 0, so we do not consider the second 
term for the covariance matrix. The quadratic estimator 
determines the tensor-to-scalar ratio using the relation 



T/S = [F-']q 

q 

F 



T/S^O ' 



- -TrfCn'C 



T/S) 



-Tr (Co '^Ct/sCq ^^t/s) ; 



(3) 



here F is the 1x1 Fisher matrix and Co can be any pos- 
itive definite hermitian matrix. In our case the Fisher 
matrix is reduced to only 1 element, since we have only 1 
parameter (T/S). For any choice of Cq the estimator will 
be unbiased, however the quadratic estimator method 
has optimum efficiency if Co is the covariance matrix it- 
self. In our particular case Cq contains the instrument 
noise which is uncorrelated and also contains the contri- 
bution of the scalar modes to the power spectrum. As we 
are interested in estimating A (^), we want to evaluate 
the Fisher matrix element F. In this case Ct/s con- 
tains just the i3-modes contribution to the power spec- 
trum. Equation ^ requires the signal and the noise to 
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Masks 




FIG. 1: The masks used for cutting the sky based on 
HEALPix resolution 3 galactic dust maps in Aitoff projection. 
Various levels of gray show the way the masks were built by 
adding areas with higher dust temperature. The units are of 
black body temperature. 
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FIG. 2: wCf ^, wCf ^ and the lensing + detector noise level 
after applying the Gaussian window for T/S ~ 1/2. wCf"^ is 
represented by the dashed line, wC^^ by the continuous line 
and the lensing + detector noise level by the thick horizontal 
line. These are the values used for resolution 6, where I My = 
196. 



be Gaussian. The gravitational wave signal is Gaussian 
because it comes from the fluctuations of a linear quan- 
tum oscillator. Instrument noise is typically Gaussian 
due to the central limit theorem; it is assumed here for 
simplicity that it is white but this assumption may be 
inappropriate for some experiments. 

The Fisher matrix error A (^) is the la standard devi- 
ation of T/S' assuming that T/S — 0. The caveats of this 
method arc that the error distribution of A (^) is non- 
Gaussian and it depends on T/ S (in particular if T/S is 
increased, the error goes up). During our calculations, 
we used the galactic dust temperature map to create 
the masks for our patches. For all the resolutions con- 
sidered in our calculations, we used the same dust map, 
sampled at a lower resolution. As we included a higher 
portion of the sky in our calculations, we kept adding 
regions contaminated more and more by the foregrounds 
in our map. They are based on masks, some of which 
are shown in Fig^ For one of the calculations, in order 
to point out the effect of having a shorter length for the 
boundary of the cut, we used cuts parallel to the galactic 
plane. 



III. NUMERICAL METHOD 

We have developed two different numerical methods 
for evaluating Eq. 13 both of them having as input the 
scalar and tensor modes generated by CMBFAST |23| 
with the RECFAST recombination routines [2ll[23|. We 
are using maps of various resolutions in the HEALPix 
|23Ll23 | pixelization. Given the resolution r, the number 
of pixels a HEALPix map has is = 12(4''). The pixel 
area is A = Att/N, and the side of the pixel is 



The Nyquist mode will be 



12n 
2~d 



(5) 



We are using both the scalar and the tensor modes with 
2 < I < l^y. l^y, and we apply a Gaussian window on 
them of the form 



w = exp 



-Trl{l + 1) 



3aN 



(6) 



here a is a Gaussian window factor. To avoid high-/ E- 
modes leaking down to low / through pixelization aliasing 
effects, we require that wCf^ ^ (jnoise ^ which setting 
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wCf^) 

LMy ^ 



10 gives 



12 In 



'Ny 



In (10) 



(7) 



and in FigI21 we can see the window function satisfying 
the above requirements applied to the E and i3-modcs. 

The exact method uses the convolution algorithm of 
Ref. to generate the Co matrix starting from Cf^ , 
(jnoise implement Eq. © for the 

Ct/s matrix 

from ■ After generating the matrices, we want 
to compute ^Tr [( 

C-^Ct/s)^]- We use the Cholesky 
method for computing the inverse of Cq and then the 
back substitution method for the multiplication with 
Cx/s- The final step is taking the square of the final 
matrix, which is done classically. As we want to evaluate 
the same quantities for different cuts of the sky, we can 
save processor time, by doing the computation progres- 
sively increasing the area which is covered. The computer 
memory is another limiting factor for our computation. 
Given that N is the number of pixels, this method needs 

(4) 

to store in the memory 3A^^/4 double precision numbers. 
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FIG. 3: wC; for two different values of tlie optical depth 
to the reionization r. The dashed line corresponds to a value 
of r = 0.07 and the continuous line to r = 0.17. The thick 
horizontal line is the lensing + detector noise level. 



to the leading order of N. With 4GB of memory avail- 
able, r = 5 is the highest resolution at which we can 
perform the exact calculation. The exact method is an 
A'''^ process, which makes the processor time required for 
higher resolutions to increase rapidly. 

In addition to the method described above, we are us- 
ing a Monte Carlo method to evaluate the Fisher matrix. 
It evaluates 



IV. RESULTS 

We have computed the Fisher matrix using Eq. 
and from this, we can calculate 



(9) 



Fr 



T/S,T/S 



(8) 



We used both methods for resolutions 4 and 5, while 
for resolution 6, we used only the Monte Carlo method. 
In order to measure more accurately the tensor-to-scalar 
ratio, a scientific mission would need to cover a higher 
fraction of the sky. We make plots showing the uncer- 
tainty in the tensor-to-scalar ratio (A (^) ) versus the sky 
fraction we need to cover in order to achieve this uncer- 
tainty. As the Monte Carlo method does not output the 
exact value of the Fisher matrix, but it only estimates it, 
every time we used the Monte Carlo method, the points 
are presented with the corresponding la error bars. The 
area of the sky used in our calculations is chosen based 
on the masks which were presented in a previous section 
(Fig©. 

For the entire sky, the uncertainty in the tensor-to- 
scalar ratio is about 1.5 x 10^^. According to the usual 
scaling A (^) cx f^f}J^, for a 70% coverage of the sky 
we should obtain A (^) = 1.8 x 10"^; instead, from our 
calculations we get A(^) =3.2x10^^, nearly a factor of 
2 degradation over the ideahzed case. For fsky > 0.7 the 
scaling is f~^^y, as it is depicted in Fig^by the straight 
diagonal line. While the scaling is less steep below that 
we still find it is steeper than expected. To achieve a 3 
a detection oiT/S = 10~^ we need 15% sky coverage, 
compared to 0.3% from idealized scaling, while for T/ S = 
10^^ the required sky coverage is 1%. 



where x is a randomly generated vector and C,^ x is 
computed iteratively, without having to take the in- 
verse of the Co matrix. We evaluate C(^^x using 
the un-preconditioned conjugate gradient method with 
the spherical harmonic transform routines described in 
Ref. [23 . We compute the mean for different realizations 
of X and use the standard deviation to obtain la errors. 
This method is more efficient than the exact approach 
for resolutions higher than 4 as it scales as (x N^/"^ (the 
number of conjugate-gradient iterations needed to solve 
C^^yi and the the number of realizations of x needed to 
compute the trace are roughly 400 and 25, respectively, 
and do not depend significantly on the resolution, but 
depend more on the number of pixels chosen within the 
mask and also on the choice of the mask). We chose a 
value of 5.0 /iK arcmin for the noise, which is the contri- 
bution of the lensing noise level. The results can be scaled 
down to lower lensing -I- detector noise if lens cleaning is 
used [5I l2^ l27l |. e.g. a 2.5 /j,K arcmin experiment would 
have an uncertainty A(T/S) that is reduced by a factor 
of 4. 
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FIG. 4: A as a function of sky fraction. The dashed 
line is computed based on the exact method for resolution 
4, while the points with error bars are based on the Monte 
Carlo method for resolution 6 The continuous lines represent 
the theoretical scaling of A The lower line corresponds 

'T\ _ r-1/2 



to the idealized scaling A ( g-j cx /^j.^ , while the upper line is 
A (X f~i^y and fits better the actual results for f^ky > 0.7. 



To interpret these results the most relevant question is 
what constraints do these results place on the required 
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level of foreground subtraction. Typically in CMB we 
remove some part of the sky where foregrounds are just 
too difficult to deal with, while in the remaining area 
some foreground subtraction method is applied to reduce 
the contamination below the required level. This level of 
course depends not only on the foregrounds, but also on 
the expected signal. To get a lower A (^) we have to map 
a larger fraction of the sky, but this means that we have 
to get into more contaminated regions of the sky, which 
requires better methods for cleaning the foreground. The 
cleaning methods are based on expected frequency scal- 
ing of dust radiation, which differs from that of CMB. 
However, frequency scaling may not be the same at ev- 
ery direction because of a different composition of dust 
grains. This leads to a decorrelation of dust maps at 
differing frequencies. This provides a fundamental limi- 
tation to foreground cleaning that is very difficult to cir- 
cumvent. Since we do not have data about the polariza- 
tion of dust radiation we will just assume different levels 
for it. We also do not know how decorrelated dust maps 
are as a function of frequency separation. Rough expec- 
tation is that dust is polarized at a few percent level, 
so we expect that cleaning at 1% should be relatively 
straightforward, at 0.1% significantly more challenging 
and at 0.01% may be impossible. The tensor-to-scalar 
ratio for which the polarized signal is equal to the po- 
larized fraction of the dust thermal emission is defined 
as: 

where p is the polarized fraction of the dust thermal emis- 
sion. We will assume we work at 90 GHz, where dust sig- 
nal is relatively low and we ignore the contribution from 
polarized synchrotron emission, which is expected to be 
similar to or slightly less than the dust polarization at 
90 GHz . Because of this many of the upcoming po- 
larization experiments are expected to operate at higher 
frequencies where dust signal is stronger relative to CMB, 
in which case our results may be overly optimistic and 
need to be rescaled appropriately. Note that our anal- 
ysis is simplistic in the sense that we only use Eq. (|10|) 
to asses the detectability prospects and do not include a 
more detailed information on how the two signals vary 
with scale. It should nevertheless give a good ballpark 
estimate of the required levels of decontamination. 

For example, if we want to detect a tensor-to-scalar 
ratio of 10^**, we would need to map 70% of the sky and 
the polarized fraction of the dust thermal emission needs 
to be less than 0.01%., while if we want to detect a tensor- 
to-scalar ratio of 10^'^, we need to map 15% of the sky 
and the polarized fraction of the uncleancd dust thermal 
emission in this region cannot exceed 0.1% (Fig[SJ). 

Besides parameters that depend on the instrument, 
there are also cosmological parameters that can infiu- 
ence the detection of i?-modes. As it was discussed in 
the Introduction section of the paper, a cosmological pa- 
rameter which has a large influence on the results of our 
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FIG. 5: The lines represent the fraction of the sky for which 
the polarized fraction of the dust thermal emission is less than 
the polarized signal, as defined in equation 1101 assuming the 
polarized fraction of the dust thermal emission to be 1% (top), 
0.1% (middle) and 0.01% (bottom). Dashed area is excluded 
if we require a detection of tensor-to-scalar ratio at a 3cr 
confidence. For example, to detect T/S — 10~^ at 3a we 
need to observe 15% of the sky and if that is chosen to be the 
cleanest region of the sky then the dust polarization in this 
region needs to be cleaned at 0.1% level of its intensity. 
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FIG. 6; A computed based on the exact method for 
resolution 4. The dashed line corresponds to an optical depth 
T — 0.07, while the continuous line corresponds to r = 0.17. 



calculations is the optical depth to the reionization, r. 
For example, if we vary the value of t from 0.17, to 0.07, 
for HEALPix resolution 4 the value of A (^) increases 
by half an order of magnitude, wile the effect of varying 
the optical depth on the i3-modes can be seen in Fig|31 
We show A (■^) for different values of optical depth r in 
FigEl 

As space missions require a large budget, the scientific 
community is also considering ground-based or balloon- 
borne experiments. They can achieve a higher angu- 
lar resolution, but their sky coverage has to suffer. We 
are also considering this case, by doing a calculation at 
HEALPix resolution 7, confined to declinations 6 < —30° 
around the terrestrial south pole. For the purpose of this 
calculation, we ignored the contribution of the reioniza- 
toin peak to the tensor modes, by considering only C^^ 
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FIG. 7: A (-^^ for HEALPix resolution 7 for a map confined 
to Earth's southern hemisphere, dechnations S < —30°. The 
graph is obtained by ignoring the reionization peak for the 
B-modes and using the Monte Carlo method, with la error 
bars shown on the graph. The continuous line is a cubic spline 
through the considered points. 

with / > 20. The results of this calculation can be seen in 
FigEI Thus, in order to be able to detect A (|) = 10"^ 
at 3cr, we need to look at 1% of the sky and from EatTUI 
we find that the polarized dust signal must be cleaned 
down to 1% of its intensity. 

V. DISCUSSION AND CONCLUSIONS 

The main conclusion we get from these results is that 
the Fisher matrix element ^t/s,t/s does not depend lin- 
early on the sky percentage considered in the map. In- 
stead, it depends more like Ft/s,t/s f%y This is 
in some way expected as the harmonic decomposition of 
the sky is nonlocal and cutting off some parts of the sky 
affects the rest of it also especially for the large scale 
where the i3-modes are more important. It is not an ex- 
act power law dependence, but as it can be seen from 
FigEI the power law could be a good approximation. We 
fitted the points from 10% to 100% in sky fraction to 
a power law. There is a disagreement between this de- 
pendence and the one from FIG. 8 of Ref. flQl. where 
Ft/s,t/s fljjy. The method proposed in Ref. [13 uses 
projected E and -B-modes, which triggers loss of infor- 
mation in the low multipoles. In the pixel space, this 
loss increases for small sky patches. To give an idea of 
the way low multipoles affect our analysis, in the case 
where the effect of the i3-modes in the reionization peak 
was not considered (FigEI), Ft/s,t/s ^ flky We con- 



sidered different pixelizations by changing the resolution 
of our HEALPix maps and we can see from Fig0| that 
this does not affect much the noise levels, leaving the re- 
sults unchanged. The way we choose to make the cut is 
mildly important in evaluating the results. If instead of 
using a cut based on the galaxy dust maps, we take a 
straight symmetric cut, parallel to the galactic equator, 
we get higher Ft/st/s values, but the change is modest 
(Fig®. 
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FIG. 8; The difference in A for two different types of 
cuts. The continuous line is obtained for the parallel equato- 
rial cut. Both results are obtained through the exact method 
for HEALPix resolution 4. 

Our results suggest that to achieve a detection of T/ S 
at 10"'^ level one needs to observe 15% of the whole sky 
as opposed to 0.3% naively expected. This makes the 
case for a space mission stronger if such levels oi T/S 
can be motivated theoretically. However, this is some- 
what dependent on the assumed optical depth, which we 
do not know yet with enough precision. Moreover, it is 
unclear if we can actually foreground clean the observed 
portion of the sky sufficiently well; as we have indicated 
the required levels are challenging indeed, in the case of 
T/S = 10~^ they require a factor of 10-100 reduction of 
polarized dust emission even if the cleanest par t of the 
sky is chosen, agreeing with the results from where 
the authors considered several specific situations. It is 
premature to give a precise verdict in favor of a space mis- 
sion or a ground-based /balloon-borne experiment, what 
is clear is that we need more data from the experiments 
currently running in order to put a better constraint on 
the parameters that influence the detection of i?-modes, 
and then decide for the type of experiment we have to 
build next. 
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